Josephson Oscillation and Transition to Self- Trapping for Bose-Einstein- Condensates 

in a Triple- Well Trap 
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We investigate the tunnelling dynamics of Bose-Einstein-Condensates(BECs) in a symmetric as 
well as in a tilted triple- well trap within the framework of mean-field treatment. The eigenenergies 
as the functions of the zero-point energy difference between the tilted wells show a striking entan- 
gled star structure when the atomic interaction is large. We then achieve insight into the oscillation 
solutions around the corresponding eigenstates and observe several new types of Josephson oscil- 
lations. With increasing the atomic interaction, the Josephson-type oscillation is blocked and the 
self-trapping solution emerges. The condensates are self-trapped either in one well or in two wells 
but no scaling-law is observed near transition points. In particular, we find that the transition from 
the Josephson-type oscillation to the self-trapping is accompanied with some irregular regime where 
tunnelling dynamics is dominated by chaos. The above analysis is facilitated with the help of the 
Poicare section method that visualizes the motions of BECs in a reduced phase plane. 

PACS numbers: 03.75.Kk, 03.75.Lm 
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I. INTRODUCTION 

Since the first realization of dilute degenerate atomic 
gases in 1995, a new epoch for studying the dynamical 
property of Bose-Einstein Condensates (BECs) comespj. 
For the dilute degenerate gases, essential dynamical 
property is included in the Gross-Pitaevskii equation 
(GPE)|3. The nonlinearity, originated from the inter- 
atomic interaction, is included in the equation through 
a mean field term proportional to condensate den- 
sity. Previously, several authors investigated the dy- 
namics of GPE for a double-wellpotential in a two- 
mode approximation^ Bl IE IE H • Novel features 
were found, e.g., the emergence of new nonlinear station- 
ary states and a variety of new crossing scenarios 0, 
nonzero adiabatic tunnelling probability^, Q, etc, to 
name only a few. Among these findings, nonlinear 
Josephson oscillation and self-trapping phenomenon are 
of most interest. As well be known, for single particle 
in a symmetric double- well, the tunnelling dynamics is 
determined by the tunnelling splitting of two nearly de- 
generate eigen-states and tunnelling time or quantum os- 
cillation period is inversely proportional to the energy 
splitting [lfj. When atomic interaction emerges, the tun- 
nelling between two-well is also observed, termed as non- 
linear Josephson oscillation[E IE Hi- However, in this 
case, the oscillation period sensitively depends on the 
initial state but has little relation to the difference be- 
tween the eigenenergies. More interestingly, with increas- 
ing the atomic interaction further (even it is repulsive), 
the Josephson oscillation between two wells is completely 
blocked, the BECs atoms in a symmetric double-well po- 
tential shows a highly asymmetric distribution as if most 
atoms are trapped in one This somehow counter- 

intuitive phenomenon is termed as self-trapping and has 
been observed in lab recentlv[Tl|. 

In the present paper, we extend to investigate tun- 



nelling dynamics for BECs in a triple- well system (l2Lll3| 
(schematically sketched as Fig^l, want to know how the 
nonlinear Josephson oscillation and self-trapping behave 
in this simplest multi-well system. This extension is not 
trivial because quantum tunnelling may happen between 
several wells simultaneously so that we expect that the 
tunnelling dynamics in the triple-well will show more in- 
teresting behavior. Moreover, the study of triple-well 
system will provide a bridge between the simple double- 
well and the multi-well, helping us understand the 'self- 
localized' phenomenon of BECs in the optical lattice[l4|. 

Technically, to investigate the dynamics of triple- well 
systems we resort to the Poicare section method[l5| that 
visualizes the motion of BECs in a reduced two dimen- 
sional phase plane. For our triple-well system, ignoring 
a total phase the dynamics is governed by a Hamiltonian 
with two freedoms. Its phase space is four-dimensional. 
However, the motions in a high-dimensional phase space 
( in our case, it is 4D) are difficult to trace. With us- 




FIG. 1: The schematic sketch of our model, (a) The sym- 
metric case(7 = 0); (b) The asymmetric case, (7, 0, —7) is the 
zero-point energy in each well respectively. 
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ing the Poicare section, we can investigate the motions 
of BEC in a reduced 2D phase plane. 

Our paper is organized as follows. In Sec. II we in- 
troduce our model and show the unusual structure of 
the eigen-energies. In Sec. Ill we investigate nonlin- 
ear Josephson oscillations with using Poincare section 
method and demonstrate diverse types of the oscillations 
for BECs. In Sec. IV we investigate the transition from 
the Josephson oscillation to self-trapping of BECs in one- 
well as well as in two-well and show an irregular regime 
characterized by chaos. Our discussions are extended to 
tilted triple- well system in Sec.V. Final section is our dis- 
cussions and conclusions. 



II. MODEL 

For the triple- well system and under mean field approx- 
imation, the wave function \I>(r, t) of GPE is the superpo- 
sition of three wave functions describing the condensate 
in each trap, i.e., 




FIG. 2: The eigen-energy levels for different interaction 
strength. We have set v = 1. 

and the corresponding canonical equations 



*(r, t) = oi(t)^i(r) + a 2 {t)<t>2{r) + az{t)fo{r). (1) 

Then the triple-well system is described by a dimen- 
sionless Schrodinger equation, 
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~~dt 
d9 1 

~dt 



= -v sin {9 1 ) ^fn{\J\ — n\ — n 3 , 



7 + -c (2ni - 2 (1 - Tlx - n 3 )) 
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with the Hamiltonian 



H = 



7 + c|ai| 2 







c a 2 



-l + c\a 3 \ 2 



(3) 



The total probability | ai | +|(X2| +l a 3| is conserved and 
is set to be unit, c is the mean field parameter denoting 
the atomic interaction and v is the coupling parameter, 
7 is the zero-point energy of the wells. The schematic 
sketch of the model is shown in FigQ] In our following 
discussions, we focus on the case of repulsive interaction 
between atoms, i.e., c > 0. 

With ignoring a total phase, the dynamics of the 
above three-level quantum system can be depicted by 
a classical Hamiltonian of two-degree freedom [^. Let 
us set that, m = |ai| 2 ,n 2 = |a2 1 2 , «-3 = |a3 1 2 , 6*i = 
argai — arg 02,6*3 = arga 3 — arga2, use the constraint 
ni + rii + n 3 — 1, we can get the classical Josephson 
Hamiltonian, 



H = 



7(774 - n 3 ) + (nf + n| + (1 - m - n 3 ) 2 ) 
—vy/1 — m — n 3 {y/n{cos9i + ^/n 3 ~cos9 3 ) , (4) 



dn 3 
dt 

d(h 
dt 



-v sin (6» 3 ) y/l-m- n 3 ^/ni, 

~1 + \f (2n 3 - 2 (1 - ni - n 3 )) 

v ( v /nT cos (6>i) + cos (fl 3 ) ^/nj) 
2^1 -n\-n 3 



(5b) 
(5c) 



v cos (9 3 ) y/l — m— n 3 



(5d) 



The fixed point or minimum energy point of the classi- 
cal Hamiltonian system (HI corresponds to the eigen-state 
of quantum system^, Il6j . To derive the analytical ex- 
pressions of these fixed points is difficult, however, nu- 
merically, we can readi ly o btain them with exploiting the 
Mathematica Software |L7| . We plot the eigenenergies as 
the function of the zero-point energy bias in Fig. 2, they 
show unusual entangled star structure for the strong non- 
linearity. 

For the weak interactions, the eigen-energy levels is 
very similar to linear case(c = 0). With increasing the 
nonlinearity(i.e., c=1.5), topological structure of the up- 
per level changes: two small loops emerge. When the 
interaction is stronger(i.e., c=2.5), the two loops will col- 
lide and form a star structure, and then two more loops 
will emerge at the middle level E%. For still stronger 
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interaction (i.e., c=4.5), the star structure of the up- 
per level entangles with the star structure in the lower 
level. However, we can still distinguish these levels be- 
cause they have different relative phases. In fact, levels 
labeled by E\ have relative phases (9i = tt, #3 = 7r), lev- 
els labeled by E2 have relative phases (#1 = tt, 63 = 0) 
or (6>i = 0,6*3 =■ tt), levels labeled by £"3 have relative 
phases (6 1 =0,9 3 = 0). 

The relation between the chemical potential and the 
above energy 

£ = M-^(M 4 + M 4 + |a 3 | 4 ), (6) 

where /1 denotes the chemical potential defined as 
(*|JJ|*>. 

In the above calculations and henceforth, for conve- 
nience we set the coupling parameter as unit, i.e., v = 1. 

III. THE SYMMETRIC TRIPLE-TRAP CASE, 

7 = 



The stability of the corresponding eigen-states can be 
evaluated by the eigen- values of the Jacobian of the clas- 
sical Josephson Hamiltonian Q). 



a 2 n d 2 n a 2 n a 2 n \ 

d 2 H d 2 H d 2 H d 2 H 

dn'f d9idn 1 dn^drii dO^dni 

d 2 H d 2 H d 2 H d 2 H 

dntdea dBtdOa dn 3 dd 3 ~d6%~ 

d 2 H d 2 H d 2 H d 2 H 

dn\dn 3 d0\dn 3 dni. d0 3 dn 3 / 



The eigen- values of the above Jacobian have their cor- 
respondence of the Bogoliubov excitation spectrum of 
BECs . Pure imaginary values indicates to a stable BECs 
state, whereas emergence of real values implies the un- 
stable for BECs and lead to a rapid production of the Bo- 
goliubov quasi-particles^8j. From calculating the above 
Jacobian matrix and making diagonalization we know 
that, the states corresponding to level En, and level E-^ 
are unstable, others are stable. 



Firstly we focus on the symmetric case, i.e., 7 = 0. 
The dependence of the energy levels on the interaction 
strength is exposed by Fig. 3. For the interaction strength 
is less than a critical value c\ — 1.56, the level structure is 
similar to its linear counterpart except for some positive 
shifts on the energy values. For c > ci, two more levels 
labeled as [E\ a ,E\b) emerge. Actually, they correspond 
to the star structure of upper level Ei in FigOand have 
relative phases of (9i = tt, 63 — tt). When the interaction 
strength is still stronger and exceeds the second critical 
value C2 = 4.06, the other two more energy levels labeled 
as {E2a, E2b) emerge. They correspond to the star struc- 
ture of mid-level E2 in FigEl an( l have the same relative 
phases as that of E2- 



A. Linear Josephson Oscillation Solution 

For the linear Josephson oscillation, i.e., c = 0, the sys- 
tem is analytically solvable. The solutions of (a\, 0,2,0,3) 
are 



01 = C 2 cos(— t + C 3 ) + Ci, (8a) 
a 2 = C 4 cos(-^=i + C 5 ), (8b) 

a 3 = -(C 2 cos(^=(t+'^)+C 3 ) + C 1 ). (8c) 



V2 V 




1 2 3 4 5 



FIG. 4: The evolutions of ni (heavy line), n2(thin line), 
n.3 (dashed line) and Si (heavy line), #3 (dashed line) in sym- 
metric wells for linear case(a) and weak interactions case(b), 

FIG. 3: When 7 = 0, the energy levels vary with the interac- with the same initial values (m = 0.4, 713 = 0.05, 6\ — 0, 63 = 

tion strength c/v, where we set v = 1. 0). 
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FIG. 5: The Poincare section at #3 = for c/v = 0.5 with 
different energy E. (a)E = -0.4, (b)E = -0.1, (c)E = 0.2, 
(d)£ = 0.5. 
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Where Ci are parameter determined by initial conditions 

C 2 cosC 3 + Ci =ai(0), (9a) 
C4COSC5 = a 2 (0), (9b) 
C 2 cosC 3 -Ci =o 3 (0), (9c) 

__^_C 2 sinC 3 = J ^a 2 (0), (9d) 
-^=C 4 sinC5=z^(a 1 (0)+a 3 (0)), (9e) 

and the constraint 

M0)| 2 + |a 2 (0)| 2 + |a 3 (0)| 2 = l. (10) 

From the above explicit expressions, we see that a%, a 2 
and et 3 vary with respect to time periodically. They share 
a common period that is inversely proportional to the 
coupling parameter, i.e., T = 2^/2n/v. Actually, the 
frequency is just the bias between eigen-energy levels. 
With initial conditions, the coefficients Ci will be fixed by 
using ©. In our case, C\ 7^ 0, so the period of population 
ni and n 3 is twice the period of n 2 , and compared with 
m, variable n 3 has a phase delay of half-period. The 
above analysis is confirmed by our numerical simulations 
as shown in FigQJa). 

B. Weak Interaction cases, c < ci 

The dynamics of this two-freedom system could be vi- 
sualized from Poincare sectional . We do this by solving 
the canonical equations l|5d|l numerically and then plot- 
ting 6*i and n\ at each time that 6* 2 — and 62 < 0. No- 
tice the total energy is conserved, therefore the Poincare 



FIG. 6: The poincare section at 63 = for c/v — 5 with 
different energy E. (&)E = 0.3, (b)E = 0.8, (c)E = 1.1, 
(d)E = 1.5, (e)E = 1.9, (i)E = 2.3. 

sections consists of a picture panel where each picture 
corresponds a fixed energy. 

For the linear case, all of motions share a common pe- 
riod exactly, and the Poincare section is some isolated 
points. With weak interactions, the periodicity will be 
destroyed, and the motions become periodic or quasi- 
periodic, corresponding Poincare section is plotted in 
FigEI where we see the section plane is full of stable 
islands. However, in this case the motion is similar to 
the linear case if they have the same initial values, as 
shown in FigQJa)0Ib). 

C. Strong Interaction Cases, c > C2 

When the interaction is strong, the nonlinear effect 
is dominant, accordingly the Poincare section is compli- 
cated, as shown in FigEJ where we see many chaotic 
region except for some stable islands. In the islands, 
the motions are periodic or quasi-periodic whereas in the 
chaotic region the motions are irregular. In order to grasp 
the dynamical property in this situation, we have simu- 
lated the motions for every regular islands numerically, 
and find that except for the oscillations like the linear or 
weak interaction shown in Fig[7fa) , there are also 

four types of new oscillations, as shown in FigJTfb) EJc), 
respectively. 

Fig[7{a) shows that the oscillation in well one is almost 
the same as that of well three except a phase delay of 
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FIG. 7: The evolutions of n\ (heavy line), n2(thin line), n 3 (dashed line) and 0i (heavy line), #3 (dashed line) in symmetric 
wells for strong interactions(c/« = 5), with initial values (a)same to FigUl (ni = 0.4, n 3 = 0.05, Oi — O,0 3 = 0); (b) (ni = 
0.49, n 3 = O.O12,0i = 0, 3 = 2.82); (c)(m = 0.15, n 3 = 0, 0i = 0, 3 = 0); (d)(m = 0.076, n 3 = 0.066, 0i = 2.54, 3 = 2.79); 
(e)(m = 0.96, n 3 = 0.02, 0i = 1.36, 3 = 1.34); (f)(m = 0.348, n 3 = 0.532, 0i = 0.817, 3 = 3.843). 



half-period. In order to compare with the linear or weak 
interaction case, we take the same initial value as that of 
Fig0] We see that their oscillations behavior is similar. 

In addition to the case shown in Fig[7Ja), the motions 
of BECs in triple- well can demonstrate very different be- 
havior. FigE^b) shows that almost all of BECs atoms 
oscillate with small amplitude in two adjacent wells, i.e., 
well one and well two. The phase d\ oscillates around 
and the phase 83 oscillates around it. The energy of these 
oscillations is closed to the eigen-energy of level labeled 
by E2 a , and the center they surrounded is near the fixed 
point corresponding to level i5 2a . As mentioned before, 
this fixed point are stable point. 

Fig[7Jc) shows that almost all of BECs atoms oscillate 
with large amplitude in well one and well two, and the 
relative phase 9\ is always oscillating around zero. These 
oscillations are regarded as oscillations in a reduced two- 



well trapped system. In fact, because n 3 is small and c/v 
is very large, we can regard the term 



-Hi = vy/l - n~i - n 3v ^cos (0 3 ) (11) 
as a perturbation. Using the generating function 

G = vg ( Ji, J 2 ) sin (<9 3 ) + J x d x + J 2 8 3 . (12) 

Where 

/ t T , Vy/l - Ji - J2VJ2 

5(Jl ' j2H c(l-J!-2J 2 ) • 
Then the Hamiltonian becomes 
1 



H' = ic (2 Jf + 2 J 2 Ji - 2 Jx + 2 Jf - 2 J 2 + 1) 



-iVl - J\ - ^2\AAC0S ($1) 



(13) 
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The new canonical variables have relations with the old 
canonical variables 



m = Ji, 

n s = J 2 + v cos (6*3)5 {-h, J2) , 
$a =0 1 +i-sin(0 3 ).9 (1 ' o) (Ji, Ja). 
$2 = 6 ! 2 + usin(6>3).9 (0 ' 1) (Ji, J 2 ). 



(14a) 
(14b) 

(14c) 

(14d) 



Since the action variable J2 is constant, the action- 
angle variables Ji,<I>i can be solved first from the new 
canonical equations l|14l) . Notice that Ji = n\ is the 
population in well one, $1 « 9\ is the relative phase of 
quantum state in well one and well two. The oscillations 
shown in Figdc) just like the zero-phase mode oscilla- 
tions in a two- well trapped system. 

FigEtd) and 0^e) show self-trapping of BECs in the 
middle well and self-trapping in one side well respectively. 
These motions have a close relations to the property of 
fixed points, and will be discussed in detail at the next 
section. 

FigEJf) shows the chaotic trajectory which corre- 
sponds to the chaotic region in the Poincare section. In 
this case, the population in each well shows irregular os- 
cillation with respect to time. 



IV. TRANSITION TO SELF-TRAPPING 



<n> 



<n> 



<n> 




A. Self- Trapping in One Well 



Self-trapping is caused by the nonlinear interactions. 
For symmetric two-well system, as we have known, self- 
trapping happens only when the interaction parameter 
excesses a critical value. By calculating the averaged 
population for same initial value a\(0) = 1,02(0) = 
with different interaction strength c/v, we show this 
in FiglSfa). The Hamiltonian used for calculation is 

12 v 



H 



cm 



c\a 2 \ 



It is clearly shown that, when 



c/v > 2 the averaged ri\ is no longer zero indicating the 
beginning of the self-trapping. Soon after that, with in- 
creasing the interaction the BECs will be trapped in well 
one completely. Because the transition corresponds to 
crossing over a separatrix from oscillation to liberation, 
at the transition point c/v = 2 the scaling law follows a 
logarithm function [l9L |20| . 

For our triple-well system, the high-dimensional phase 
space permits the existence of chaos and the smooth 
movement of the fixed points to the boundary ( the latter 
point will be clearly shown in Fig. 9 ). So we expect that 
the transition to self-trapping in triple-well system will 
show distinguished property from that of double- well sys- 
tem. To demonstrate it, we calculate the time averaged 
population for different interactions with the initial con- 
ditions 722(0) = 1 and ni(0) = 1, denoting initial BECs 



FIG. 8: The average of ni,ri2 for different interactions c/v 
in two- well system(a) and triple- well system(b,c), with ini- 
tial value (a)m(O) = 1, (b)ra 2 (0) = 1, (c)m(O) = 1. The 
schematic sketch of the potential is shown in the figures. 



uploaded in middle well and left-hand well, respectively. 
The results are shown in FiglH^t)) Etc) - 

For the linear case, i.e., c = 0, in FigHJb), substituting 
(ai = 0, a?, = 1, as = 0) to iJSJl, we can get one set of Ci 

1 jr 

C\ = 0, C% = — — — , C3 = — , C4 = 1, C5 = 0. 
Then from Eq.® we have the averaged populations 

<n 1 ) = (n 3 ) = |C7 2 | 2 /2=i <n 2 ) = |C7 4 | 2 /2 = ^ 

Likewise, in FiglSJc), substituting (a± = 1,02 = 0,03 = 
0) to we can get one set of Ci 

Cl = r C2 = r C3 = > C4 = -^ C5= r 

From Eq. (JSJ we have the averaged populations 

1 - (n 2 ) 3 



(n 2 ) = |C 4 |V2=-,(n 1 ) = (n3) , 



7 



1.0 



c 

TO 



LU 

uf0.5 



c 
o 



O 

a o.o 







1- 






A 

CM 

C 
+ 

V 




VjV^V 




5^a*Z n/(E 1a ) (n 2 *'(E 1a )) 


0- 
1" 









1 



c/v 



FIG. 9: The population n[ (heavy line) and n^thin line) for 
the levels labeled by Ei,Ei a and Eu in FigEl 



For the weak interaction, i.e., c < ci, the averaged 
populations are still similar to the linear case. However, 
when the interaction strength is close to the critical point 
Ci, the averaged populations changed dramatically. For 
FigEfb), the averaged population (712) increases mono- 
tonically and smoothly to unit. No scaling law are ob- 
served. This is due to the fact that in the 4D phase space 
the fixed point can move smoothly to the boundary with- 
out through any bifurcation. For FigHJc), the averaged 
populations become turbulence near the critical point, 
this is a result of chaotic trajectory in the phase space. 
Meanwhile, (ni) and (713) are no longer equal. When the 
interaction is larger than 2.25, the averaged population 
(ni) becomes smooth and tends to unit rapidly. 

For triple-well system, when the interaction strength 
exceeds the critical value c\, as mentioned before, new 
levels Ei a ,Eib will appear, accordingly new fixed points 
will emerge, and the phase space will tend to be divided 
into several subspace around the stable fixed points. In 
Fig|5]we plot populations n\ and n\ for the levels labeled 
by Ei,Ei a and En, as a function of c/v. Corresponding 
to eigenenergy E\ a there are two eigen-states, denoted 
by n* and n*' respectively. The same thing happens to 
level En,. Recall that Ei,Ei a are stable levels. When 
the interaction is very strong, the averaged population 
(71-2) of the motions of BECs uploaded initially in the 
middle well is close to n^-Ei), while (ni) of the motions 
of BECs uploaded initially in well one is close to nj (Ei a ), 
as shown in FigOJb),Fig|BIc). 



B. Self- Trapping in Two Wells 

In the above, we have investigated the self-trapping 
of BECs in single well. In this part, we will investigate 



+ 

V 





(b) 



V e 3=° 



-6,-63=0 



■e l -e 3 = It /2 



c/v 
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FIG. 10: The mean value of (a)(ni + 712) and (b)(ni + 713) 
with initial value (a)(ni = 0.5, 712 = 0.5, 713 = 0), (b)(ni = 
0.5, 712 = 0, 713 = 0.5) for different relative phase. 



whether the BECs atoms can be self-trapped only in two 
wells. 

Considering the interference, the initial relative phase 
should be very important. So we calculate the mean value 
of (711+712) as a function of c/v for different phase Q\ with 
initial value [ri\ = 0.5, 77.2 = 0.5,773 = 0), and the mean 
value of (n! + 77 3 ) as a function of c/v for different phase 
81 — 9 3 with initial value [n\ = 0.5, n 2 = 0,n 3 = 0.5), 
respectively. The main results are shown in Fig llOl 

It is shown that, with increasing the atomic interaction 
BECs will be trapped in the two wells where it is initially 
uploaded and the Josephson oscillation can be completely 
blocked. Both cases also suggest that the relative phase 
can dramatically influence the transition to self-trapping 
for BECs. In Fig llOf a). when relative phase is zero, the 
Josephson oscillation and the self-trapping can emerge 
alternately. Whereas, for case of Fig HUT b). it value of the 
relative phase gives a robust self-trapped BECs. In both 
cases, we also see the occurrence of the chaos making the 
curves look irregular . Interestingly, the onset of chaos 
can also be controlled by the relative phase, e.g., the 
vanishment of the relative phase will reduce the chaos 
and make the BECs safely turn from oscillation state to 
self-trapping state as shown in Fig llOf b). 
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V. TILTED TRIPLE- WELL, 7 / 



In this section we extend the above discussions to the 
tilted triple-well system. In this system, the diverse type 
of Josephson oscillations also emerge around the eigen- 
states. In light of the energy spectra Fig. 2, we can read- 
ily find out the parameter regime for different kinds of 
oscillations. So, we will focus on the transition to self- 
trapping that is of more interest. We want to see how 
the transition is influenced by tilting the wells. 



A. Linear Oscillation Solutions 

For the linear case, i.e., c = 0, the system is analyti- 
cally solvable. So the solutions of (01,02,03) are 



1 v 2 1 

ai = -Ca cosh/ — + 7 2 i -I- C5) -I o 

2 V V 2 ' b> 2(^+ 7 2 ) 



{lC 2 cos{\j— + 1 H + C 3 )+vC l ), (15a) 



a.2 



f +7 2 2 
7C1), 



(15b) 



1 



1 



a 3 = -Ca cos(\ V 7 2 i + C 5 ) - 

2 4 V V 2 ' b> 2(f + 7 2 ) 



("/C 2 cos(\l—+ 1 H + C 3 )+vC 1 ). (15c) 



<n> 



<n> 



<n> 




Where C, are complex integral constant and like the 
symmetric case, determined by initial conditions. 

It is clear that a\ , a 2 and 03 vary with respect to time 



periodically. The period T = 2-kU Pf + 7 2 



B. Transition to Self- Trapping 

When the interaction is strong we still observe the self- 
trapping of BECs in the tilted system. We plot the aver- 
aged populations in their dependence of the interaction 
strength with choosing parameter 7 = — 1 and initial con- 
ditions n\ = 1, n,2 = 1 and n 3 = 1 in Fig lllf al lllf b^ and 
Hir e 1 ) respectively. The schematic sketch of the wells is 
shown in Fig^fb). 

For Fig lllf a). when c/v = 0, with the initial conditions 
(di = 1, a 2 = 0, a 3 = 0), we derive one set of Cj 



d = v/2, C 2 = iv/V2, C 3 = arccos(^^), 

v 

C 4 = 



y/v 2 + 2 7 2 



G 5 = — arccos(— ). 



FIG. 11: The average of rii for different interactions c/v in 
asymmetric wells 7 = — 1 with initial value (a)m(0) = 1, 
(b)na(0) = 1, (c)n 3 (0) = 1. 



By integrating |ai(£)| 2 with respect to time and mak- 
ing average over time we obtain the averaged population 
analytically from ljl5|l . 

(m) = 17/24, (n 2 ) = 1/4, (n 3 ) = 1/24. 

With increasing the nonlinearity, the averaged popula- 
tion in well one (ni) decrease at first, and then become 
turbulence, indicating chaotic motions. When the inter- 
action parameter is larger than 3.9, the averaged popu- 
lation jumps up and tends to unit soon after. 

In Fig lllf b). initial condition is (ai = 0, a 2 = l,a 3 — 
0). Similarly the averaged populations are readily ob- 
tained for the linear case, (ni) = 1/4, (n 2 ) — 1/2, (n 3 ) — 
1/4. With increasing the nonlinearity, the averaged popu- 
lation (n 2 ) increases smoothly at beginning, passes a tur- 
bulence interval [2.3,3.6], and then jumps up and tends 
to unit. Another interesting phenomenon in the process 
is that the averaged population in well three is always 
larger than that of well one in the presence of the non- 
linearity, even though in this case the zero-point energy 
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of well three is obviously greater than that of well one. 
This somehow counterintuitive phenomenon is clearly the 
consequence of the nonlinearity. 

In Fig lllf cl. when c/v — 0, with the initial conditions 
(ai = 0, a,2 = 0,03 = 1), the averaged populations have 
some correspondence to that of Fig lllf a]) due to behind 
symmetry, i.e., (m) = 1/24, (n 2 ) = 1/4, (n 3 ) = 17/24. 
With increasing the nonlinearity, the averaged popu- 
lation in well three (n^) increases monotonically and 
smoothly to unit. 

With comparing Fig^l to FigEl we find that, the 
smooth transition of the BECs in middle well to self- 
trapping in the symmetric triple- well is broken by tilting 
the wells, and lifting the well three makes BECs smoothly 
transit to self-trapping states without losing their stabil- 
ity. Fig. 11 also shows that, BECs in higher wells are eas- 
ily self-trapped. From the above discussions, we conclude 
that the transition to self-trapping of BECs in triple-well 
systems can be effectively controlled by tilting the wells. 

VI. CONCLUSIONS 

We have presented a comprehensive analysis of the 
tunnelling dynamics for BECs in a triple-well trap both 
numerically and analytically. Diverse energy levels are 
demonstrated. Behind these unusual level structures, we 
reveal many new types of nonlinear Josephson oscillation. 



We also study the self-trapping of BECs in one-well as 
well as in the two- well and investigate the transition from 
nonlinear Josephson oscillation to self-trappings. Dis- 
tinguished from the double-well case, no scaling law is 
observed at the transition and the transition may be ac- 
companied by a irregular regime where the motions are 
dominated by chaos. We also find that the transition can 
be effectively controlled by the relative phase between 
wells and tilting the wells. In the present experiments, 
the double-well is realized in the optical traps with using 
a blue-detuned light to form a barrier. With the same 
technique, the triple-well is also possibly realized in the 
optical traps. We hope our theoretical discussion will 
stimulate the experiments in the direction. 
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